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FOREWORD 
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Inc.,  Armament  Systems  Division,  7225  Northland  Drive,  Brooklyn  Park,  Minnesota 
55428,  for  the  White  Oak  Laboratory,  Naval  Surface  Warfare  Center  (NSWC), 

10901  New  Hampshire  Avenue,  Silver  Spring,  Maryland  20903-5000,  under  contract 
N60921-86-C-0249. 

This  effort  was  conducted  during  the  period  from  November  1986  to 
September  1987.  The  authors  would  like  to  thank  Paula  Walter,  NSWC  program 
manager,  and  Drs.  Kibong  Kim  and  Frank  Zerilli  for  helpful  technical  discussions 
during  the  course  of  this  program. 
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was  performed  with  Honeywell  Independent  Development  funding. 
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SECTION  1 
INTRODUCTION 


This  report  provides  user  instructions  for  the  1987  version1  of  the  EPIC-2/ 
REZONE  coae  for  use  with  the  1986  version  of  the  EPIC-2  code.  Use  is  restricted  to 
two-dimensional  geometry  (axisymmetric,  plane  strain,  plane  stress),  and  only 
triangular  elements  can  be  rezoned. 

The  primary  motivation  for  this  work  was  the  desire  to  be  able  to  rezone  com¬ 
putations  involving  explosive-metal  interaction  and  subsequent  jet  formation.  It 
should  be  useful  for  other  applications. 

There  are  two  primary  reasons  for  rezoning.  One  is  to  improve  accuracy  and  the 
other  is  to  decrease  computing  time  through  use  of  a  larger  integration  time 
increment.  When  triangular  elements  become  very  distorted  (long  and  slender)  they 
may  not  be  able  to  properly  represent  gradients  alonj*  the  longer  dimension  of  the 
element,  and  this  can  lead  to  artificial  stiffness.  Similarly,  during  high  distortion  the 
minimum  altitude  of  the  element  can  become  very  small  and  this  decreases  the 
integration  time  increment,  which  is  bounded  by  the  sound  speed  transit  time  across 
the  element.  Rezoning  allows  the  user  to  replace  distorted  elements  with  compact 
elements;  decreasing  the  larger  dimension  gives  greater  accuracy,  and  increasing  the 
shorter  dimension  gives  a  larger  integration  time  increment. 

There  are  two  separate  functions  performed  in  the  EPIC-2/REZONE  code.  The 
first  is  to  delete  the  distorted  grid  and  to  replace  it  with  a  better  grid  (more  compact 
elements).  The  second  function  is  to  transfer  the  physical  properties  (node  and 
element  variables)  from  the  old  grid  to  the  new  grid.  This  second  function  is  per¬ 
formed  automatically  without  user  intervention. 

The  remainder  of  this  report  provides  a  brief  description  of  the  approach,  user 
instructions  and  examples. 


1  Johnson,  G.R.  and  Stryk,  R.A.,  User  Instructions  for  the  EPIC-2  Code. 
AFATL-TR-86-51,  Honeywell  Inc.,  Defense  Systems  Division,  Edina,  MN, 
Sep  1986. 


1/2 


NSWCTR  88-382 


SECTION  2 
APPROACH 


This  section  provides  a  brief  overview  of  how  the  EPIC-2/REZONE  code  works. 
The  following  is  a  list  of  the  primary  functions  performed  by  the  code. 

•  The  EPIC-2/REZONE  code  reads  the  restart  file  which  must  have  been 
previously  generated  by  the  EPIC-2  code. 

•  The  element  variables  from  the  restart  file  are  transferred  to  nodal 
quantities.  For  example,  the  nodal  pressure  is  the  average  of  the  element 
pressures  for  all  elements  attached  to  a  specified  node. 

•  Portions  of  the  old  grid  (from  restart  file)  are  deleted  and  the  new  grid  is 
generated.  The  node  and  element  variables  from  the  old  grid  are  saved  in 
temporary  arrays. 

•  Each  new  node  searches  through  the  old  grid  to  find  the  appropriate  old 
element  which  contains  the  position  of  the  new  node. 

•  The  velocities  and  accelerations  for  the  new  node  are  linearly  interpolated 
from  the  old  nodal  velocities  and  accelerations  of  the  three  old  nodes  which 
define  the  old  element  which  contains  the  new  node. 

•  The  center  position  of  each  new  element  searches  through  the  old  grid  to 
find  the  appropriate  old  element  which  contains  the  center  position  of  the 
new  element. 

•  The  variables  for  the  new  element  are  linearly  interpolated  from  the  old 
nodal  quantities  of  the  old  element  variables. 

•  The  new  grid,  with  node  and  element  variables  interpolated  from  the  old 
grid,  is  written  on  a  restart  file.  This  new  restart  file  is  then  available  to  be 
read  by  the  EPIC-2  code  for  further  computations. 

It  should  be  noted  that  there  is  no  absolute  conservation  of  mass,  momentum  or 
energy.  If,  however,  the  outline  of  the  new  grid  closely  coincides  with  the  outline  of 
the  old  grid,  then  these  system  parameters  should  be  closely  matched  (within  a  few 
percent).  Significant  discrepancies  in  the  outlines  of  the  old  and  new  grids  can  lead  to 
significant  discrepancies  in  the  mass,  momentum  and  energy.  The  user  should 
compare  the  summaries  of  the  system  data  output  for  the  grid  before  and  after 
rezoning,  to  ascertain  the  accuracy  of  rezoning. 
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SECTION  3 

USER  INSTRUCTIONS 


This  section  provides  user  instruction  for  the  EPIC-2/REZONE  code.  The 
function  of  this  is  to  read  a  restart  file,  define  a  new  geometry,  transfer  physical 
properties  to  the  new  geometry,  and  write  a  restart  file.  If  the  rezoning  performed 
with  this  program  requires  changes  in  slide  lines,  etc.,  these  changes  must  be  made 
with  the  Change  Drop  and  Change  Add  options  in  EPIC-2. 


INPUT  DATA  FOR  EPIC-2/REZONE 

The  input  data  necessary  to  perform  the  rezoning  are  shown  in  Figures  1 
through  3  and  Reference  1.  The  descriptions  which  follow  are  for  the  data  in 
Figures  1  through  3.  The  units  used  must  be  the  same  as  those  used  for  the 
EPIC-2  preprocessor  input. 

It  is  possible  to  intellect  user  comments  into  the  data  by  use  of  a  $  character.  If 
the  $  is  in  the  first  column  of  the  card,  that  entire  card  is  ignored  as  input  data.  If  the 
$  is  beyond  the  first  column  in  the  card,  then  the  $  and  all  data  to  the  right  of  the  $ 
are  ignored. 

Before  continuing,  note  the  following  definitions: 

•  A  "deleted  node”  is  a  node  from  the  original  grid  (before  rezoning)  which 
has  been  deleted.  There  are  some  instances  when  the  coordinates  of  a 
deleted  node  are  used. 

•  An  "old  node”  is  a  node  from  the  original  grid  (before  rezoning)  which  has 
not  been  deleted  and  remains  a  part  of  the  grid  after  rezoning  occurs. 

•  A  "new  node”  is  a  node  which  is  introduced  into  the  rezoned  grid. 

•  A  "previously  generated  new  node”  is  a  node  which  was  previously 
generated  during  the  current  rezoning  activity.  There  are  some  instances 
when  a  "new  node”  may  make  use  of  the  coordinates  of  a  "previously 
generated  new  node,”  when  two  adjacent  shapes  are  being  generated. 

Description  Card  (A80) 

TITLE  =  Title  for  report. 
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Miscellaneous  Card  (715) 

CYCLE  =  Cycle  number  on  restart  file  which  will  be  rezoned. 

NPROJ  =  0  will  allow  rezoning  of  projectile  only. 

1  will  allow  rezoning  of  target  only. 


DESCRIPTION  CARO  (A 80) 

(""TITLE  - - 


MISCELLANEOUS  CARO  (715) 


INDIVIDUAL  DELETED  HOPES  CARO  -  FOR  NON  ■  0  (1615) 


NODE  GROUP  CARDS -FOR  NGN  >0(315) 

I  H'A  |  NNO  |  wc 

MDMOUAL  DELETED  ELEMENTS  CARO  -  FOR  NGL  -  0  (1515) 

I  Li  |  u  i ....  i  u.  kxxxxxVxxxVVVVVVVVVVVV^^^ 

ELEMENT  GROUP  CARD  -  FOR  NGL  >  0  (315) 

|T.a-|  iM  |  wt 

NODE  DATA  CARDS  AND  ADDITIONAL  NOOE  DATA  CAROS  •  AS  REQUIRED 
(DESCRIBED  IN  FIGURE  2  OF  THIS  REPORT  AND  IN  FIGURE  3  OF  AFATL-TR-86-51) 

|  BLANK  CARD  ?  ENOS  NOOE  DATA 

ELEMENT  DATA  CARDS  ANO  ADDITIONAL  NOOE  DATA  CAROS  •  AS  REQUIRED 
(DESCRIBED  IN  FIGURE  3  OF  THIS  REPORT  ANO  IN  FIGURE  4  OF  AFATL-TR-55-51) 

|  BLANK  CARO  }  ENOS  ELEMENT  DATA 

MERGE  CARO  (FI 0.0) 

LTOLMRQ  ^\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\^ 


FIGURE  1.  EPIC-2/REZONE  INPUT  DATA 
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F0UR-8OED  SHAPE  NOOE  CARO  (515.5XF10.0) 


CORNER  NOOE  OSHNITON  CAROS  •  4  REQUIRE  (215  ;2F  10.0) 


SDE  OPTIONS  CARD  (415) 


SIDE  NOOE  COORDINATE  CAROS  •  AS  NECESSARY  (Sf  10.0) 


w.e.fr 

CORNER  NOOE  DEFINITION  CARO  (SEE  ABOVE) 

SIDE  NOOE  COORDINATE  CAROS  •  AS  NECESSARY  (SEE  ABOVE) 
OLD  SIOE  NOOES  CARO(S)  •  AS  NECESSARY  (SEE  ABOVE) 
GENERATED  siDE  NOOES  CARO(S)  •  AS  NECESSARY  (SEE  ABOVE) 


FIGURE  2.  ADDITIONAL  NODE  INPUT  DATA 
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FOUR-StOED  SHAPE  S. EM  ENT  DESCRIPTION  CARD  (615) 
I~r  I  N1  |MAfl  \—~\  NcoqTiiAY 


THREE-810  ED  SHAPE  ELB4BI?  DESCRIPTION  CARO  (SIS) 
|  »  |  ni  |maTL  |  lT 


TRANSITION  B.  EM  ENTS  DESCRIPTION  CARO  (1615) 

|  *  1  N1  | MAH  |  LI  |  N2  |  N3  |  N4  |  NS  |  N6  |  N7  |  N8  |  NO  |  N10  |  Nil  |  N12  |  N13  | 


FIGURE  3.  ADDITIONAL  ELEMENT  INPUT  DATA 


Note:  It  is  possible  to  rezone  nodes  and/or  elements  in  both  the  projectile  and 
the  target  in  a  single  rezone  run.  If  it  is  necessary  to  do  this,  then  two 
subsequent  rezone  runs  must  be  made. 

NGN  =  Number  of  groups  of  elements  to  be  deleted.  If  NGN  =  0  then  deleted 
nodes  are  read  individually. 

NGL  =  Number  of  groups  of  elements  to  be  deleted.  If  NGL  =  0  then  deleted 
elements  are  read  individually. 

NNDEL  =  Number  of  nodes  to  be  deleted. 

NLDEL  =  Number  of  elements  to  be  deleted.  Only  elements  with  the  same 
material  number  can  be  deleted. 

PRINT  =  0  will  provide  minimal  printing. 

1  will  print  new  nodes  and  elements. 

2  will  provide  extensive  printing  for  program  checkout. 

Individual  Deleted  Nodes  Card  (1615) 

The  deleted  nodes  are  input  individually  when  NGN  =  0. 

N1...NN  =  Nodes  to  be  deleted. 

Node  Group  Card  (315) 

Use  one  card  for  each  group  of  nodes  to  be  deleted. 

NIG  =  First  node  in  group. 

NNG  =  Last  node  in  group.  May  be  left  blank  if  group  has  only  one  node. 

INC  =  Increment  between  nodes.  May  be  left  blank  if  increment  is  1 . 
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Individual  Deleted  Elements  Card  (1615) 

The  deleted  elements  are  input  individually  when  NCL  =  0. 

L1...LN  =  Elements  to  be  deleted. 

Element  Group  Card  (315) 

Use  one  card  for  each  group  of  elements  to  be  deleted. 

L1G  =  First  element  in  group. 

LNG  =  Last  element  in  group.  May  be  left  blank  if  group  has  only  one  element 
INC  =  Increment  between  elements.  May  be  left  blank  if  increment  is  1. 
Node  Data  Cards 


The  node  data  cards  for  EPIC-2,  as  described  in  Reference  1,  are  acceptable  as 
well  as  the  cards  described  in  Figure  2.  (The  Additional  Node  Data  Cards  of  Figure  2 
will  be  defined  later).  There  is  no  Scale/Shift/Rotate  card  (as  used  in  Reference  1)  so 
all  coordinates  given  must  be  the  actual  coordinate.  A  blank  card  ends  the  node  data 
input. 

Element  Data  Cards 


The  element  data  cards  for  EPIC-2,  as  described  in  Reference  1 ,  are  acceptable 
as  well  as  the  cards  described  in  Figure  3.  (The  Additional  Element  Data  Cards  of 
Figure  3  will  be  described  later.)  A  blank  card  ends  the  element  data  input. 

Merge  Card  (FI 0.0) 

The  node  merging  function  provides  a  way  of  combining  adjacent  shapes  to 
cover  a  complex  region.  Each  shape  can  be  generated  using  the  node  numbering 
required  by  the  node  generator.  Adjacent  shapes  may  require  different  node  numbers 
for  the  same  node  location.  This  requirement  can  be  met  by  generating  a  node  for 
each  node  number  needed.  After  the  node  and  element  generators  have  been  used, 
then  the  multiple  nodes  (at  a  single  location)  are  reduced  to  single  nodes  by  the  node 
merging  function. 

TOLMRG  =  Greater  than  zero  to  merge  nodes  which  are  close  together.  All  nodes 
in  a  box  (TOLMRG  on  each  side)  are  combined  to  become  the  same 
node.  The  node  number  entered  first  in  the  node  input  cards  will  be 
used  for  the  combined  node.  Only  the  newly  entered  nodes  are 
merged.  TOLMRG  =  0  gives  no  node  merging. 

Additional  Node  Data  Cards 


The  node  data  cards  described  in  Reference  1  can  be  used  to  describe  the  new 
grid  but  three  new  generators  have  been  added  to  take  advantage  of  the  availability 
of  the  old  grid.  The  new  generators  can  use  the  coordinates  of  deleted  nodes  for  the 

Eosition  oi  a  new  node.  The  generators  can  also  find  the  edge  of  the  deleted  grid 
etween  two  deleted  nodes  and  use  this  line  for  the  side  of  a  new  shape.  When  this 
edge  option  is  used,  only  the  approximate  position  of  the  deleted  end  nodes  need  to  be 
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given.  With  these  features  it  is  possible  to  lay  out  a  new  grid  using  only  a  geometry 
state  plot  (See  Appendix  A) . 


Four-Sided  Shape  Node  Card  (515.  5X.  FI 0.0) 

The  four-sided  shape  is  shown  in  Figure  4.  Six  cards  are  required  for  each  four¬ 
sided  shape  generated.  Additional  cards  are  needed  for  certain  options.  The  corner 
nodes  are  denned  in  a  counterclockwise  direction  with  the  first  corner  node  being  Nl. 
The  sides  are  defined  in  a  counterclockwise  direction.  The  first  side  starts  at  the  first 
corner  node  and  ends  at  the  second  corner  node.  The  fourth  side  starts  at  the  fourth 
corner  node  and  ends  at  the  first  corner  node.  The  nodes  are  numbered  in  layers.  The 
first  layer  starts  at  the  first  corner  node  and  goes  across  to  the  fourth  corner  node. 

The  following  layers  go  into  the  interior  of  the  shape  with  the  last  layer  going  from 
the  second  to  the  third  corner  nodes.  This  shape  can  be  made  continuous  with 
adjacent  shapes  by  aligning  border  nodes  and  using  the  merge  option.  The  RELAX 
option  (described  later)  provides  a  method  of  more  evenly  spacing  the  interior 
primary  nodes  of  this  shape.  The  node  placement  algorithm  gives  good  results  for  a 
rectangular  shape  but  questionable  placement  for  distorted  shapes,  especially  those 
with  curved  sides.  When  questionable  node  placement  occurs,  the  RELAX  option  can 
improve  the  placement.  The  border  nodes  are  not  moved.  The  interior  primary  nodes 
are  moved  to  be  at  the  average  position  of  the  other  four  connected  primary  nodes. 
The  secondary  nodes  are  generated  later. 


7 

Nl 

NC 

NL 

IRFIX 

RELAX 


Identification  number  for  four-sided  node  shape  geometry. 

Number  of  the  first  node  of  the  four-sided  shape.  This  should  be  selected 
so  the  generated  new  nodes  do  not  have  the  same  node  numbers  as 
remaining  old  nodes. 

Number  of  nodes  along  the  second  and  fourth  side.  Does  not  include  the 
secondary  (crossed  triangle)  nodes. 

Number  of  nodes  along  the  first  and  third  side.  Does  not  include  the 
secondary  (crossed  triangle)  nodes. 

Radial  restraint  option.  If  IRFIX  =  1,  all  nodes  on  the  axis  of  symmetry 
are  restrained  radially.  Leave  blank  for  no  restraint. 

Less  than  or  equal  to  zero  for  no  node  motion.  Greater  than  zero  to  move 
the  interior  primary  nodes  to  the  average  position  of  the  adjacent  four 
primary  nodes.  The  primary  nodes  are  moved  in  several  sweeps.  Each 
sweep  takes  the  interior  nodes  in  numeric  order,  adjusts  the  node 
position,  and  sums  the  node  movements.  The  sweeps  terminate  when 
the  total  node  movement  for  the  last  sweep  is  less  than  RELAX.  The 
relaxed  four-sided  shape  will  have  more  evenly  spaced  nodes.  The 
secondary  nodes  are  placed  at  the  average  position  of  the  four  primary 
nodes  in  the  composite  element  (four  crossed  triangles). 


Corner  Node  Definition  Cards  (215,  FI 0.0) 


Four  cards  are  required,  one  for  each  corner  node.  The  first  corner  node  is  Nl. 
The  other  three  corner  nodes  are  defined  in  counterclockwise  order. 
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ELEMBfT  [LI  I 


FIGURE  4.  FOUR-SIDED  SHAPE  GEOMETRY 


NOLD 

NGEN 

RCNR 


The  old/deleted  node  to  use  for  coordinates.  Leave  blank  if  you  do  not 
want  this  option. 

The  previously  generated  new  node  to  use  for  coordinates.  The 
duplication  of  node  definitions  can  be  removed  by  the  merge  step. 
Leave  blank  if  you  do  not  want  this  option. 

R  coordinate  of  corner  node. 


ZCNR  =  Z  coordinate  of  corner  node. 


Note:  The  R  and  Z  coordinates  are  ignored  when  either  the  NOLD  or 
NGEN  option  is  used. 


Side  Options  Card  (415) 


There  are  five  options  (0,1, 2, 3, 4)  for  each  of  the  four  sides  (ISl,  IS2,  IS3,  IS4). 


151 

152 

153 

154 


=  0  Will  equally  space  the  side  nodes  on  a  straight  line  between  the 
corner  nodes. 

=  1  Will  equally  space  the  side  nodes  on  the  edge  of  the  deleted  region.  If 
the  corner  nodes  are  not  already  defined  as  old/deleted  nodes  on  the 
edge,  then  the  corner  node  (whose  coordinates  were  input  as  RCNR 
and  ZCNR)  is  moved  to  the  closest  old/deleted  edge  node. 

=  2  Will  position  the  side  nodes  at  coordinates  given  on  the  following 
Side  Node  Coordinate  Cards. 
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=  3  Will  place  the  side  nodes  at  the  coordinates  of  the  old/deleted  nodes 
given  on  the  following  Old  Side  Nodes  Card(s). 

=  4  Will  place  the  side  nodes  at  the  same  coordinates  as  the  previously 
generated  new  nodes  given  on  the  following  Generated  Side  Nodes 
Card(s). 

Note:  For  side  options  2, 3,  and  4,  an  additional  card  (or  cards)  is 
required  to  define  side  coordinates  and/or  side  node  numbers.  The  side  1 
option  should  be  input  first  (if  needed),  followed  by  sides  2, 3,  and  4. 

Side  Node  Coordinate  Cards  (8F10.0) 

These  cards  must  be  input  in  counterclockwise  order. 

R1...RN  =  R  coordinates  of  the  border  nodes  on  the  side.  (Do  not  include  corner 
nodes.)  One  or  more  cards  as  necessary. 

Z1...ZN  =  Z  coordinates  of  the  border  nodes  on  the  side.  (Do  not  include  corner 

nodes.)  One  or  more  cards  as  necessary. 

Old  Side  Nodes  Cards  (1615) 


These  cards  must  be  input  in  counterclockwise  order. 

N01...N0N  =  Node  number  of  old/deleted  nodes  where  side  nodes  are  to  be  placed. 

(Do  not  include  corner  nodes.)  One  or  more  cards  as  necessary. 

Generated  side  Nodes  Cards  (1615) 

These  cards  must  be  input  in  counterclockwise  order. 

NG1...NGN  =  Node  numbers  of  previously  generated  new  nodes  where  the  side 
nodes  are  to  be  placed.  (Do  not  include  corner  nodes.)  One  or  more 
cards  as  necessary. 

Three-Sided  Shape  Node  Cards  (315. 5X.  15. 5X,  FI 0.0) 

The  three-sided  shape  is  shown  in  Figure  5.  Five  cards  are  required  for  each 
three-sided  shape  generated.  Additional  cards  are  needed  for  certain  options.  The 
corner  nodes  are  defined  in  a  counterclockwise  direction  with  the  first  corner  node 
being  Nl.  The  sides  are  defined  in  a  counterclockwise  direction.  The  first  side  starts 
at  the  first  corner  node  and  ends  at  the  second  corner  node.  The  third  side  starts  at 
the  third  corner  node  and  ends  at  the  first  corner  node.  The  shape  is  meant  to  be 
roughly  triangular.  The  first  and  third  sides  have  the  same  number  of  nodes, 

NRING.  The  second  side  has  more  nodes,  2*  (NRING  -1)  -1.  The  total  number  of 
nodes  is  2*  (NRING  *  NRING  -1).  The  shape  is  composed  of  two  regions.  The  first 
region  starts  at  the  first  corner  and  extends  to  the  last  node  before  the  second  and 
third  corners.  The  remaining  region,  a  strip  along  the  second  side,  is  a  strip  of  crossed 
triangles.  The  node  numbering  is  also  done  in  two  regions.  The  first  region  is 
numbered  as  described  for  the  four-sided  shape.  The  crossed-triangle  strip  along  the 
second  side  is  numbered  in  a  direction  generally  from  corner  three  towards  corner 
two.  The  secondary  nodes  are  numbered  first  and  the  outside  primary  nodes  are 
numbered  last.  The  highest  numbered  node  is  corner  two. 
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FIGURE  5.  THREE-SIDED  SHAPE  GEOMETRY 


®  —  Identification  number  for  three-sided  shape  geometry. 

N1  =  Number  of  the  first  node  of  the  three-sided  shape. 

NRING  =  Number  of  nodes  on  two  shorter  sides  of  shape  (Sides  1  and  3).  Must  be  3 
or  larger. 

Radial  restraint  option.  If  IRFIX  —  1,  all  nodes  on  the  axis  of  symmetry 
are  restrained  radially.  Leave  blank  for  no  restraint. 

RELAX  =  Less  than  or  equal  to  zero  for  nodemotion.  Greater  than  zero  to  move 
the  interior  primary  nodes  to  the  average  position  of  the  adjacent  four 
primary  nodes.  The  primary  nodes  are  moved  in  several  sweeps.  Each 
sweep  takes  the  interior  nodes  in  numeric  order,  adjusts  the  node 
position,  and  sums  the  node  movements.  The  sweeps  terminate  when 
the  total  node  movement  for  the  last  sweep  was  less  than  RELAX.  The 
relaxed  three-sided  shape  will  have  more  evenly  space  nodes. 

Corner  Node  Definition  Cards  (215.  E10.0) 

Three  cards  are  required,  one  for  each  corner  node.  The  first  corner  node  is  Nl. 

The  other  two  nodes  are  defined  in  counterclockwise  order. 

NOLD  =  The  old/deleted  node  to  use  for  coordinates.  Leave  blank  if  you  do  not 
want  this  option. 

NGEN  =  The  previously  generated  new  node  to  use  for  coordinates.  The 

duplication  of  node  definitions  can  be  removed  by  the  merge  step. 

Leave  blank  if  you  do  not  want  this  option. 
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RCNR  =  R  coordinate  of  corner  node. 

ZCNR  =  Z  coordinate  of  corner  node. 

Note:  The  R  and  Z  coordinates  are  not  needed  when  either  the  NOLD  or 
NGEN  option  is  used. 

Side  Option  Card  (315) 

The  choice  of  five  options  (0,1  >2,3,4)  for  each  of  the  three  sides  (IS1,  IS2,  IS3)  is 
the  same  given  for  the  four-sided  shape. 

Side  Node  Coordinate  Cards 


Same  as  given  for  the  four-sided  shape. 
Old  Side  Node  Coordinate  Cards 


Same  as  given  for  the  four-sided  shape. 
Generated  Side  Node  Coordinate  Cards 


Same  as  given  for  the  four-sided  shape. 

Regrid  Line  of  Nodes  Card  (415, 2X.  311,15) 

Two  cards  are  required  if  a  single  node  is  to  be  generated  and  three  or  more 
cards  are  required  for  more  than  one  node.  The  nodes  may  be  numbered  consecu¬ 
tively  or  incremented  by  INC.  New  nodes  can  be  generated  using  the  coordinates  of 
old/deleted  nodes. 


9  =  Identification  number  for  regrid  line  of  nodes  geometry. 

N1  =  Number  of  the  first  node. 


NN 

INC 

IR 


Number  of  the  last  node  of  the  line.  Leave  blank  if  a  single  node  is  to  be 
generated. 

Node  number  increment  between  corresponding  nodes.  Leave  blank  for 
consecutive  numbering  or  if  a  single  node  is  to  be  generated. 

Radial  restraint.  If  IR  =  1,  all  nodes  will  be  restrained  in  the  radial 
direction.  Leave  blank  for  no  restraint. 


IZ  =  Axial  restraint.  IfIZ  =  1,  all  nodes  will  be  restrained  in  the  axial 

direction.  Leave  blank  for  no  restraint. 


IT  =  Theta  restraint.  If  IT  =  1,  all  nodes  will  be  restrained  in  the  0  direction. 

Leave  blank  for  no  restraint. 

IS1  =  Side  option  same  as  described  for  the  Side  Option  Card  above.  Used  only 
when  more  than  2  nodes  are  generated. 
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Corner  Node  Definition  Card 

Same  as  for  the  three-  and  four-sided  shapes.  One  corner  card  required  for 
single  node  generation  and  two  corner  cards  required  for  two  or  more  nodes 
generated. 

Side  Node  Coordinate  Cards 

Same  as  the  three-  and  four-sided  shapes.  Required  only  when  generating  three 
or  more  nodes  and  ISl  =  2. 

Old  Side  Nodes  Coordinate  Cards 

Same  as  for  the  three-  and  four-sided  shapes.  Required  only  when  generating 
three  or  more  nodes  and  ISl  =  3. 

Generated  Side  Nodes  Coordinate  Cards 

Same  as  for  the  three-  and  four-sided  shapes.  Required  only  when  generating 
three  or  more  nodes  and  ISl  =  4. 


ADDITIONAL  ELEMENT  INPUT  DATA 

The  element  data  must  be  input  in  a  manner  consistent  with  the  Node  Data 
Cards  of  Reference  1  and/or  the  Additional  Node  Data  Cards  described  in  this  report. 
In  addition  to  the  four-sided  and  three-sided  shapes,  a  Transition  Element  shape  is 
included  to  assist  in  gridding  transition  regions.  A  description  of  the  additional 
element  input  data,  for  the  three  new  shapes,  is  given  in  Figure  3. 

Four-Sided  Shape  Element  Description  Card  (615) 

One  card  is  required  for  each  four-sided  shape,  as  shown  in  Figure  4.  The  ele¬ 
ments  are  numbered  consecutively,  beginning  with  the  layer  along  side  4  (from 
corner  1  to  corner  4),  and  ending  with  the  layer  along  side  2  (from  corner  2  to 
corner  3).  The  number  of  elements  generated  is  4*NCOL*NLAY. 

7  =  Identification  number  for  four-sided  shape  geometry. 

N1  =  Number  of  first  node  in  four-sided  shape. 

MATL  =  Material  number  for  the  elements.  Must  be  identical  to  the  deleted 
material. 

Ll  =  Number  of  the  first  element  in  the  four-sided  shape. 

NCOL  =  Number  of  columns  of  composite  elements  (four  crossed  triangles)  along 
sides  2  and  4. 

NLAY  =  Number  of  layers  of  composite  elements  (four  crossed  triangles)  along 
sides  1  ind  3. 
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Three-Sided  Shape  Element  Description  Card  (515) 

One  card  is  required  for  each  three-sided  shape,  as  shown  in  Figure  5.  The 
element  numbering  follows  the  pattern  of  node  numbering.  The  first  region  is  num¬ 
bered  first  and  the  strip  along  the  second  side  is  numbered  from  the  third  corner 
towards  the  second  corner.  The  number  of  elements  generated  is  4  *  (NRCOL  * 
NRCOL-1). 

8  =  Identification  number  for  three-sided  shape  geometry. 

N1  =  Number  of  first  node  in  three-sided  shape. 

MATL  =  Material  number  for  the  elements.  Must  be  identical  to  the  deleted 
material. 

LI  =  Number  of  first  element  in  three-sided  shape. 

NCOL  =  Number  of  composite  elements  (four  crossed  triangles)  along  sides  one 
and  three. 

Transition  Elements  Description  Card  (1615) 

One  card  is  required  for  each  group  of  transition  elements,  as  shown  in 
Figure  6.  The  group  has  a  central  node  and  form  3  to  12  elements  around  the 
central  node. 

9  =  Identification  number  for  transition  geometry. 

N1  =  Number  of  the  central  node. 

MATL  =  Material  number  for  the  elements.  Must  be  identical  to  the  deleted 
material. 


N3 


FIGURE  6.  TRANSITION  SHAPE  GEOMETRY 
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LI  =  Number  of  the  first  element.  The  first  element  will  be  composed  of  the 

nodes  Nl,  N2,  and  N3.  The  second  element  will  be  Nl,  N3,  and  N4.  The 
last  element  will  be  Nl,  NK,  and  N2  where  NK  is  the  last  non-blank 
perimeter  node  given. 

N2..N13  =  Numbers  of  the  perimeter  nodes.  At  least  3  nodes  must  be  given.  The 
order  of  the  nodes  must  be  counterclockwise.  Newly  generated  nodes 
and  old  (undeleted)  nodes  may  be  used. 


PROGRAM  STRUCTURE  AND  FILE  DESIGNATION 

A  hierarchy  chart  for  EPIC-2/REZONE  is  shown  in  Figure  7.  Many  of  the 
subroutines  are  identical  to  those  of  EPIC-2.  Table  1  contains  a  list  of  the  modified 
and  unmodified  EPIC-2  subroutines. 

The  file  designations  are  as  follows: 

IN  =  5  Input  file 

IOUT  =  6  Output  file 

ITAPIN  =  9  Restart  tape  read  by  REZONE 

ITAPOT  =  10  Restart  tape  generated  by  REZONE 

TABLE  1 .  EPIC-2  SUBROUTINES  USED  IN  EPIC-2/REZONE 


Modified 

Unmodified 

ECORE 

ECOUNT 

ESHAPE 

NEWMAS 

NBCORE 

NFIX 

NCORE 

RECALL 

NSHAPE 

SAVE 

SDATA 

INSTRUCTIONS  FOR  CHANGING  PROGRAM  DIMENSIONS 

All  the  arrays  which  have  their  dimension  changed  are  in  COMMON  state¬ 
ments  with  the  dimension  given  by  a  parameter.  Each  COMMON  statement  and  its 
associated  PARAMETER  statement  is  in  a  VAX  style  INCLUDE  file.  A  complete 
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FIGURE  7.  HIERARCHY  CHART  FOR  EPIC-2/REZONE 
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list  of  all  the  COMMON  statements  can  be  found  in  the  BLOCK  DATA  subprogram 
RZINIT.  A  program  dimension  can  be  changed  by  editing  the  PARAMETER  state¬ 
ment  in  an  include  file  and  recompiling  the  program.  A  comment  explaining  the 
meaning  of  each  parameter  follows  eacn  PARAMETER  statement. 

All  instructions  for  changing  the  program  dimensions  of  EPIC-2  applyhere  also 
except  the  DATA  statement  in  subroutine  GEOM  has  been  moved  to  RZINIT.  The 
arrays  in  EPIC-2/REZONE  must  be  large  enough  to  contain  the  problem  handled  by 
EPIC-2.  A  safe  way  to  approach  this  is  to  have  corresponding  arrays  in  both 
programs  the  same  size. 


EXAMPLE  PROBLEMS 

Figure  8  shows  three  computations  of  a  copper  cylinder  impacting  a  rigid 
surface.  In  Case  A  there  is  no  rezoning,  in  Case  B  there  is  a  rezone  at  lOps,  and  in 
Case  C  there  are  rezones  at  lOps  and  20ps.  In  all  three  cases,  the  computations  were 
performed  for  a  total  of  80ps. 
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RIGID  SURFACE 


CASE-A 

•  NO  REZONING 

•  8842  CYCLES 


*s  2. 


CASES 

•  REZONE  AT  IOjiS 

•  3796  CYCLES 


CAS.E.C 

•  REZONE  AT  10(iS  AND  20^S 

•  2040  CYCLES 


FIGURE  8.  CYLINDER  IMPACT  EXAMPLE  PROBLEM 
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It  can  be  seen  that  the  deformed  lengths,  the  maximum  deformed  diameters,  the 
strain  contours  and  maximum  strains  are  nearly  identical  for  all  three  cases.  The 
most  significant  difference  occurs  for  the  number  of  integration  cycles.  Here,  the 
number  of  cycles  is  greatly  decreased  for  the  rezone  runs,  and  this  significantly 
decreases  the  required  computing  time. 

EPIC-2  Preprocessor  innut  data  for  the  cylinder  impact  example  are  shown 
in  Figure  9  and  EPIC-2/REZONE  input  data  for  the  rezone  at  lOps  are  shown  in 
Figure  10.  The  intent  of  this  example  is  to  demonstrate  the  conservation  of 
properties  after  rezoning  and  the  reduction  in  computing  time. 

An  example  of  an  explosive-metal  interaction  computation,  with  subsequent  jet 
formation,  is  shown  in  Figure  11.  For  this  problem,  the  severe  distortions  make  it 
desirable  to  change  the  nodal  connectivity  such  that  there  are  more  elements  along 
the  length  of  the  jet.  The  resulting  rezone  then  increases  both  the  computational 
accuracy  and  the  integration  time  increment. 

EPIC-2  Preprocessor  input  data  for  thejet  formation  example  are  shown 
in  Figure  12,  ana  the  corresponding  EPIC-2/REZONE  input  data  are  shown  in 
Figure  13.  Some  features  of  the  grid,  before  and  after  rezoning,  are  shown  in 
Figure  14. 
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101  2CYLINDER  IMPACT  EXAMPLE  PROBLEM 
10  1  0  0  0  0  1 
1010  999.  COPPER 


1.0 

1.0 

2  1 

5  32 

0  1 

125984 

0. 

.125984 

1.0 

1.0 

2  1 

1  4 

31  5 

0. 

-8937. 

0. 

0.5E-7 

0. 

1. 

0  0 

1. 

0.1E-8 

0  0 

0  0 

000005 

1111. 

0.0 

000010 

1111. 

0.0 

000015 

1111. 

0.0 

000020 

1111. 

0.0 

000040 

1111. 

0.0 

000060 

1111. 

0.0 

1. 

0. 

0.0 

1.27559  0.  1. 

0. 


1 


0.9  .000060  1.0 

1 

1 

1 

1 

1 

1 

1 


FIGURE  9.  EPIC-2  INPUT  DATA  FOR  CYLINDER  IMPACT  EXAMPLE  PROBLEM 


CYLINDER  REZONE  AT  10  MICROSECONDS 
$CYCL  NPRJ  NGN  NGL  NNDL  NLDL  PRNT 


209 

0 

1 

1 

284  496 

0 

$  NIG 

NLG 

INC 

1 

284 

1 

$  L1G 

LLG 

INC 

1 

496 

1 

$  FOUR-SIDED 

1  SHAPE  NODES 

S  7 

Nl 

NC 

NL 

IRFX///// 

RELAX 

7 

1 

5 

32 

1 

0.001 

$NOLD 

NGEN 

RCNR 

ZCNR 

1 

280 

284 

5 

$  IS1 

IS2 

IS3 

IS4 

0 

1 

1 

1 

$  BLANK  CARD  (NOTE  COLUMN  1  TO  $  IS  BLANK)  FOR  END  OF  NODES 
$  FOUR-SIDED  SHAPE  ELEMENTS 
$  7  N1  MATL  LI  NCOL  NLAY 

7  1  1  1  4  31 

$  BLANK  CARD  (NOTE  COLUMN  1  TO  $  IS  BLANK)  FOR  END  OF  ELEMENTS 
$  TOLMRG 
0.0 


FIGURE  10.  EPIC-2/REZONE  INPUT  DATA  FOR  CYLINDER 
IMPACT  EXAMPLE  PROBLEM 
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FIGURE  11.  JET  FORMATION  EXAMPLE  PROBLEM 
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1  2  EXAMPLE  JET  FORMATION  PROBLEM 


1  0 

1  0 

0 

0 

1 

1  0 

1  0 

1.0 

999. 

41  0 

1.0 

1.0 

2  1 

21  4 
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1.0 
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1.0 

2  201 

21  18 

4 

1 
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2.0 
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0.0 
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0.0 

0.0 

0.0 
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0.0 

0.0 

0.0 
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2 
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1 
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-1 
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1 
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1111.0 

0.0 

1 

0.51E-4 

1111.0 

0.0 

1 

0.6E-4 

1111.0 

0.0 

1 

0.7E-4 

1111.0 

0.0 

1 

0.1E-0 

0.0 

0.0 

1 

2.30 

2.45 


2.45 

0.0 


0.0 


FIGURE  12.  EPIC-2  INPUT  DATA  FOR  JET  FORMATION  EXAMPLE  PROBLEM 
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JET  FORMATION  EXAMPLE  PROBLEM 


SCYCL 

NPRJ 

NGN  NGL  NNDL 

NLDL  PRNT 

788 

0 

1  1  144 

240  0 

$  NIG 

NLG 

INC 

1 

144 

$  L1G 

LLG 

INC 

1 

240 

S  REGRID  LINE  OF  NODES 

9  9 

N1 

NN  INC  IRZT 

IS1 

9 

1 

8  1  100 

0 

$NOLD  NGEN 

RCNR 

ZCNR 

124 

S  THREE-SIDED  SHAPE  NODES 

FOR  REGION  A 

$  8 

N1 

NRNG/////I RF I X/////  RELAX 

8 

7 

5  1 

0.01 

SNOLD  NGEN 

RCNR 

ZCNR 

0 

3 

12 

1 

$  IS1 

IS2 

IS3 

0 

1 

0 

S  POUR-SIDED  SHAPE  NODES  fOR  REGION  B 


9  7 

Nl 

NC 

NL 

iRrx///// 

RELAX 

7 

45 

5 

5 

1 

0.01 

SNOLD 

NGEN 

3 

RCNR 

ZCNR 

16 

5 

12 

$  IS1 

IS2 

IS3 

IS4 

0 

0 

1 

0 

$  FOUR-SIDED 

SHAPE  NODES  FOR  REGION 

c 

S  7 

Nl 

NC 

NL 

IRFX///// 

RELAX 

7 

104 

5 

3 

1  0 

0.01 

SNOLD 

NGEN 

e 

RCNR 

ZCNR 

124 

9 

138 

16 

$  IS1 

1S2 

IS3 

IS4 

0 

1 

0 

0 

S  FOUR-SIDED 

SHAPE  NODES  fOR  REGION 

D 

S  7 

Nl 

NC 

NL 

IRrx///// 

RELAX 

7 

127 

3 

3 

0 

0.00 

SNOLD 

NGEN 

RCNR 

ZCNR 

16 

138 

141 

18 

$  IS1 

IS2 

IS3 

IS4 

0 

1 

0 

1 

S  REGRID  LINE  Of  NODES  IN  REGION  E 

S  9 

Nl 

NN 

INC 

IRZT  IS1 

9 

140 

0 

0 

0  0 

SNOLD 

NGEN 

RCNR 

ZCNR 

80 

$  FOUR-SIDED 

SHAPE  NODES  FOR  REGION 

r 

S  7 

Nl 

NC 

NL 

IRrx///// 

RELAX 

7 

141 

3 

2 

0 

0.01 

SNOLD 

NGEN 

RCNR 

ZCNR 

19 

142 

144 

21 

$  IS1 

IS2 

IS3 

IS4 

FIGURE  13.  EPIC-2/REZONE  INPUT  DATA  FOR  JET  FORMATION 
EXAMPLE  PROBLEM 
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0  3  0  3 

9  OLD  MODES 
X43 
20 

<M0TB  COtUHH  1  TO  $  IS  BLANK >  TOR  END  or  MODES 
I  THREE-SIDED  SHARE  ELEMENTS  rOR  REGION  A 
9  t  Ml  NATL  LI  MCOL 

•  7  1  1  4 

$  FOUR— SIDED  SHARE  ELEMENTS  FOR  REGION  B 
9  7  Ml  HATL  LI  MCOL  NLAY 

9  45  2  <1  4  4 

9  POOR-SIDED  SHAVE  ELEMENTS  FOR  REGION  C 
9  7  Ml  HATL  LI  MCOL  NLAY 

7  104  3  157  4  2 

9  FOUR-SIDED  SHARE  ELEMENTS  rOR  REGION  D 
9  7  Ml  HATL  LI  MCOL  NLAY 

7  127  4  189  2  2 

9  TRANSITION  ELEMENTS  FOR  REGION  E 

9  «  NATe  -lie  *3  M4  1,5  N6  N7  **9  NJ 

9  140  5  205  129  134  139  146  141 

9  POOR-SIDED  SHARE  ELEMENTS  FOR  REGION  F 
9  7  Nl  HATL  LI  MCOL  NLAY 

7  141  6  210  2  1 

3  SARD  (N0IE  C0LUMN  1  TO  $  IS  BLANK)  FOR  END  OF  ELEMENTS 

9  TOLHRG 

0.01 


FIGURE  13.  (Cont.) 
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FIGURE  14.  GRID  IN  THE  JET  BEFORE  AND  AFTER  REZONING 
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SECTION  4 

CONCLUSIONS  AND  RECOMMENDATIONS 


User  Instructions  have  been  provided  for  the  EPIC-2/REZONE  code.  The 
primary  emphasis  has  been  for  jet  formation  computations,  where  it  is  anticipated 
that  improvements  can  be  realized  for  both  accuracy  and  computing  time. 
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APPENDIX  A 

EPIC-2  GEOMETRY  STATE  PLOTS  WITH  NODE  NUMBERS 

The  EPIC-2  Postprocessor  for  State  Plots  has  been  expanded  to  print  node  numbers  on  the  ge¬ 
ometry  plots.  This  added  capability  is  very  helpful  when  rezoning.  The  modified  data  input  for 
geometry  plots  are  shown  in  Figure  A-l.  The  variable  PRINT  has  been  added. 

PRINT  »  0  for  no  node  numbers. 

■  1  to  print  node  numbers  next  to  the  nodes  specified  by  the  NODE  option. 

The  drawn  node  numbers  use.  the  same  size  characters  as  those  in  the  title.  The  user  must 
choose  the  size  of  the  graph  to  give  readability.  The  node  point  is  used  as  the  lower  left  corner  of 
the  box  in  which  the  number  is  printed.  Figure  A-2  shows  an  example  of  this  option,  as  used  for 
the  jet  formation  example  problem 

Figure  A-3  gives  the  modifications  to  EPIC-2/POST  1  in  the  VAX  EDIT  format,  which  will  al¬ 
low  the  node  numbers  to  be  printed.  The  variable  IPRINT  which  is  read  in  the  POST  1  main 
program  must  be  passed  to  subroutine  GEOM  and  EXGEOM.  In  GEOM  the  calls  to  subroutine 
MAPRZ  to  plot  projectile  and  target  nodes  are  replaced  by  calls  to  a  new  subroutine  MAPNOD. 
The  subroutine  MAPNOD  is  listed  in  Figure  A-4. 


FIGURE  A-1.  MODIFIED  INPUT  DATA  FOR  THE  GEOMETRY  PLOT  CARD 
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INSERT  3|C  fOSTl  raviaad  foe  noda  nuabar  plotting  87/08/19 

REPLACE  81 | C  *  SPRINT  -  1  HILL  PRINT  NODAL  QUANTITIES  ON  CONTOUR  PLOTS 

INSERT  62)C  *  AND  PLOT  NODE  NUMBERS  NEXT  TO  NODES  POR  IPNT.NE.O 

INSERT  )C  *  ON  GEOMETRY  PLOTS 

SUBSTITUTE/I PNT/IPNT> IPRINT/  330 

SUBSTITUTE/I PNT/IPNT, IPRINT/  100384 

SUBSTITUTE/I PNT , TITLE/Z PNT ,  I  PR  I NT/  2902421 

8UB8TITUTE/RMAX/  TITLE, RMAX/  2902422 

REPLACE  2902423 ;C  REV  87/OS/ll 

SUBSTITUTE/I PNT/I PNT, IPRINT/  2902486 

SUBSTITUTE/I PNT/I PNT, IPRINT/  2902707 

REPLACE  2902709 ;C  REV  87/05/11 

INSERT  2902727|  INTEGER  IPRINT 

INSERT  | C  *  IPRINT  FLAGS  NODE  PRINTING  OPTION 

SUBSTITUTE/  PLOT  THE  NODE//  2902766 
INSERT  2902767  )C  PLOT  THE  NODE 

INSERT  ) C  IP  IPRINT  SET  THEN  PLOT  NODE  NUMBER 

SUBSTITUTE/  PLOT  THE  NODE//  2902768 
INSERT  2902769 )C  PLOT  THE  NODE 

INSERT  )C  IF  IPRINT  SET  THEN  PLOT  NODE  NUMBER 

REPLACE  2902999)  CALL  MA?NOD( NNX , LINTYP , NSYM , IPRINT ) 

REPLACE  2903010;  CALL  MAPNOD(NNX, LINTYP, NSYM, IPRINT) 

EXIT  POSTl . POR/SEQUENCE 


FIGURE  A-3.  MODIFICATIONS  TO  EPIC-2/POST  1  TO  ALLOW  FOR  NODE 
NUMBER  PLOTTING 
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SUBROUTINE  MAPNOD (NA.LINTYP , NSYM , X PRINT ) 

•  REVISED  87/05/11 

•  NAl'S  A  NODE  ONTO  THE  GRAPH  AND  OPTIONALLY  PRINT  NODE  NUMBER 
INTEGER  NA,LXNTYP, NSYM, SPRINT 

•  CALLS  WINDOW 

•  CALLED  BY  GEOM 

•  LATEST  REVISION  KAY  1907 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

•IP  PLOT  TO  RXQBT  OP  Z  AXIS  THEN 
•IP  PLOT  NOT  ROTATED  THEN 
PLOT  NORHAL  NODE 
•IP  SPRINT  SET  THEN 

PRINT  NODE  NUMBER  WHEN  ON  GRAPH 
•ENDXP 

•ELSE  AXIBS  ROTATED 
PLOT  ROTATED  RODE 
•ZP  SPRINT  SET  THEN 

PRINT  NODE  NUMBER  WHEN  ON  GRAPH 
•BNDZP 

•ENDXP  ROTATED 
•BNDZP  RIGHT 

•IP  PLOT  TO  LSPT  OP  Z  AXIS  THEN 
•SP  PLOT  MOT  ROTATED  THEN 
PLOT  HIRROR  NODE 
•ZP  SPRINT  SET  THEN 

PRINT  NODE  NUMBER  WHEN  ON  GRAPH 
•ENDZP 

•ELSE  AXZES  ROTATED 
PLOT  ROTATED  HIRROR  NODE 
•ZP  SPRINT  SET  THEN 

PRINT  NODE  NUMBER  WHEN  ON  GRAPH 
•ENDXP 

•ENDZP  ROTATED 
•ENDXP  LBPT 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

LOGICAL  AXROT 

•  AXROT  INDICATED  AXIAL  ROTATATXON  <R,Z  TO  2,-R) 

LOGICAL  LEPT 

•  LEPT  INDICATES  TO  PLOT  TO  LEPT  OP  Z  AXIS  (NEGATIVE  R) 

REAL  R 

•  R  IS  THE  RADIAL  COORDINATE  OP  A  NODE 
LOGICAL  RIGHT 

•  RIGHT  INDICATES  TO  PLOT  TO  RIGHT  OP  Z  AXIS  (POSITIVE  R) 

REAL  X 

•  X  CONTAINS  HORIZONTAL  INFORMATION 

•  X( 3 ) ,X( 4 )  ARE  SCALEXNG  INFO  PROM  SCALE  CALL  EARLIER 
REAL  Y 

•  Y  CONTAINS  VERTICAL  INFORMATION 
REAL  Z 

•  Z  IS  THE  Z  COORDINATE  OP  A  NODE 

HXMD  IS  THE  MAXIMUM  NUMBER  FOR  NODES 
PARAMETER  MXND-1000 

COMMON /MODE/R ( HXMD > , Z  ( MXND ) ,  ROOT ( MXND ) ,  ZDOT ( MXND ) ,  TOOT ( MXND ) , 

C  AHASS(MXND)  ,PHAS3(MXND)  , COUNT ( MXND )  , PNEW( MXND)  , NODE ( MXND ) , 

C  SPAREN  ( MXND  ),NODP(  MXND)  ,  GRX(MXND)  ,R0(MXND)  ,  VERT  (MXND) 

COMMON  /CHAP/  AXROT, LEPT, RIGHT, X< 4 ), Y( 4 ), HLONG,VLONG, 

S  RHZN,HMAX,VMXN,VMAX 
C 

RMA  -  MA 

C  *IP  PLOT  TO  RIGHT  OP  Z  AXIS  THEN 
XP( .NOT. RIGHT) GOTO  1 
C  *IP  PLOT  NOT  ROTATED  THEN 

I P ( AXROT ) GOTO  2 
C  PLOT  NORMAL  NODE 


FIGURE  A-4.  LISTING  OF  SUBROUTINE  MAPNOD  TO  ALLOW  FOR  NODE 
NUMBER  PLOTTING 
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CALL  WXNDOW(R(NA)  ,Z(NA)  ,R(NA)  , Z(NA)  , LINTYP , NSYM ) 

C  *IF  XRRXNT  SET  THEN 

XF( I PRINT. EQ.l) THEN 

XF(R(MA) .GE.HMIN  .AMD.  R(NA) .LE.BMAX  .AMD. 

1  2(MA) .GE.VHXH  .AND.  Z ( NA ) . LE . VMAX ) THEN 

XPAGE  -  (R(MA)-X( 3 ) )/X( 4 ) 

YPAGE  -  (Z(MA)-Y(3) )/Y( 4 ) 

CALL  MUHSERI  XPACE ,  YPAGE  ,0.14,  RMA ,  0 . 0 , -1 ) 

END  IP 
END  IP 
GOTO  3 

•ELSE  AXXES  ROTATED 
PLOT  ROTATEO  NODE  (Z(NA) ,-R(NA) ) 

2  CALL  WXNDOW(ZINA), -RINA), Z(NA), -RINA), LINTYP, NSYM) 
•IF  X PRINT  SET  THEN 

XP( XPRXMT. EQ.l )THEN 

XFI Z (NA)  .OE.HMXN  .AND.  Z ( NA)  . LE. HMAX  .AND. 

X  -R(NA) .GE.VKXN  .AND.  -R(NA) .LE. VMAX) THEN 

XPAGE  -  ( Z ( NA ) -X ( 3 ) ) /X ( 4 ) 

YPAGE  -  |-R(NA)-Y(3>)/Y<4) 

CALL  NUMBER! XPAGE, YPAGE, 0.1 4, RNA, 0. 0,-1) 

END  IF 
ENDXF 

•ENDXP  ROTATED 

3  CONTINUE 
•ENDXP  RIGHT 

•IP  PLOT  TO  LEFT  OP  Z  AXIS  THEN 
1  XP( .NOT. LEFT) GOTO  11 

•XP  PLOT  NOT  ROTATED  THEN 
X  P { AXROT ) GOTO  12 
PLOT  MIRROR  NODE  (-R(NA) . 3(NA) ) 

CALL  WINDOW ( — R ( NA )  ,  Z ( NA )  , -R ( NA )  ,  Z  ( NA )  .LINTYP  , NSYM) 
•IP  IPRINT  SET  THEN 
*F( X PRINT. EQ. 1 )THEN 

IF(-R(NA) . GE.HMIN  .AND.  -R(NA) .LE.HMAX  .AND. 

1  Z(NA) .GE.VHXN  .AND.  Z ( NA) . LE. VMAX )THEN 

XPAGE  -  (-R(NA)-XI 3 ) )/X( 4 ) 

YPAGE  -  ( Z ( NA) -Y( 3 ) )/Y( 4 ) 

CALL  NUMBER( XPAGE , YPAGE ,0.14, RNA, 0 . 0 , -1 ) 
ENDXP 
ENDIP 
GOTO  13 

•ELSE  AXXES  ROTATED 

PLOT  ROTATED  MIRROR  NODE  ( Z ( NA) , R( NA)  ) 

12  CALL  WINDOW(Z(NA),R(NA),Z(NA), RINA), LINTYP, NSYM) 

C  *IP  IPRINT  SET  THEN 

I F 1 1  PRINT .EQ.l) THEN 

IP(ZINA) .GE.HMIN  .AND.  Z( NA) . LE. HMAX  .AND. 

1  RINA) .GE.VMXN  .AND.  R( NA) . LE. VMAX) THEN 

XPAGE  -  (Z|NA)-X(!) )/X(4) 

YPAGE  -  ( R( NA )-Y( 3 ) )/Y( 4 ) 

CALL  NUMBER! XPAGE , YPAGE ,0.14, RNA, 0 .0,-1) 
ENDIP 
ENDIP 

C  *  END I P  ROTATED 

13  CONTINUE 

C  * ENDIP  LEFT 

11  RETURN 
END 


FIGURE  A-4.  (Corn.) 
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